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^ ■ Abstract 

We present the status of the development of a dedicated high density ultra-cold neutron 
(UCN) source dedicated to the gravitational spectrometer GRANIT. The source employs 
superthermal conversion of cold neutrons to UCN in superfluid helium. Tests have shown 
! that UCN produced inside the liquid can be extracted into vacuum. Furthermore a ded- 

icated neutron selection channel was tested to maintain high initial density and extract 
only neutrons with a vertical velocity component v± < 20 cm/s for the spectrometer. 
£T) . This new source would have a phase-space density of r uCN « 0.18 cm -3 (m/s) -3 for the 

spectrometer. 

od : 

o : 
o . 

Introduction 

> ■ 

•i-H . 

The solutions for Schrodinger's equation for a neutron bouncing on a reflecting horizon- 
tal surface in the earth gravitational field are given by Airy functions [lj. This textbook 
example of bound energy states in a linear potential has been demonstrated experimen- 
tally at the high flux reactor of the Institut Laue Langevin [El [3] . A new gravitational 
spectrometer, GRANIT [5j, is being built to investigate these quantum states further, and 
to induce resonant transitions between gravitationally bound quantum states. Applica- 
tions are a refined measurement of the electrical charge of the neutron, the search for 
the axion [4J, and other additional forces beyond the standard model. Experiments at the 
present UCN source are limited by counting statistics and systematic effect^]. In these 
conference proceedings the advances of a dedicated UCN source, based on superthermal 
conversion of cold neutrons in superfluid helium are described. 



"Corresponding author: Tel.: +33 4 76 20 70 27; fax: +33 4 76 20 77 77 E-mail address: schmidt-w@ill.fr 
1 For a discussion of systematic effects concerning these measurements see G. Pignol (these proceedings). 
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General aspects of Helium-4 based superthermal 
converters 



The dominant process in the conversion is the excitation of a single phonon inside the 
superfluid helium. Cold neutrons with a wavelength around 8.9 A (k = 27r/A* = 0.7 A -1 ), 
i.e. 1.0 meV kinetic energy, can be scattered down to the ultra-cold energy range by 
emission of one single phonon [6] . Multiphonon processes may also occur, depending on 
the incoming neutron spectrum pj, [8]. The resulting saturated UCN density, 

PUCN =P'T, (1) 

is determined by production rate density P and neutron storage time r inside the converter 
volume. The rate 

1111 

l/r=- + + — + , (2) 

Tf3 7"wall 7"up 7"abs 

includes contributions from all existing loss channels (neutron beta decay, losses due to 
wall collisions, up-scattering due to phonons, nuclear absorption by 3 He- impurities in the 
helium). As the absorption cross section of 4 He is zero there is no absorption inside a 
pure converter. 

The equation of detailed balance, 

a (E VCN -^E*) = -f- exp (- E *~f c A ■ a (E* — > £ UCN ) , (3) 

relates the cross sections of a two level system, as our one-phonon excitation process, de- 
pendent on temperature. For T — > K the up-scattering cross section will be arbitrarily 
small. In pure 4 He only two contributions to the total storage time remain; the neutron 
beta decay and losses due to wall collisions. The UCN production rate density is defined 
as the conversion rate of neutrons to energies below the Fermi potential of the walls of 
the converter volume and is thus dependent on the wall material. The production rate 
density expected for the 1-phonon process in a beryllium (Vp = 252 neV) coated converter 
volume (V F = 252 neV) is Pi = (4.55 ±0.25) • lO^d^/dA^s^cm- 3 , with the differential 
flux at A* = 8.9 A given in cm" 2 s _1 A _1 [9j. 

Source concept 

GRANIT and its source will be positioned on the neutron beam H172 on level C of 
the high flux reactor at the ILL (see fig.Q]). A crystal monochromator, positioned 12 m 
downstream from the cold source, will reflect 8.9 A neutrons under a take-off angle of 
9 = 61.2°. Between monochromator and converter volume a 3 m long neutron guide 
with m — 2 supermirror coating will be installed. The reduction in background will 
outweight by far the decrease in conversion rate due to the imperfect monochromator 
and the omission of multi-phonon processes. A nitrogen cooled beryllium filter will be 
installed further upstream of the monochromator to reduce high energy neutrons at the 
monochromator position. As source we will employ an existing apparatus [10] modified 
for our new requirements; a continuous flux of UCN within a narrow phase space element. 

A sketch of the converter volume with heat screens and extraction is shown in figure[2j 
The UCN will diffuse from the converter volume through a highly specular guide into 
an intermediate volume. This second volume at room temperature allows a vertical tube 
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Figure 1: Setting of the monochromator, the source, and GRANIT inside level C 
of the ILL. (For more details see: www.ill.fr/pages/science/imgs/PlanInstILL.gif) 

extraction from the converter. Thus avoiding cryogenic difficulties of a direct horizontal 
extraction with foils which has posed problems in the past [TT] and reduces background 
from directly scattered cold neutrons. Furthermore it evenly distributes all neutrons to 
the semidiffuse channel [12J. Which is used to select neutrons with small vertical energy 
E < mgh < 20 peV for the spectrometer whereas all other neutrons are reflected back 
or are lost on wall collisions. This method will increase the storage time inside the 
intermediate volume and the converter volume, as more than 85% of the neutrons incident 
on the channel entrance will be reflected back. Thus the UCN density is increased which 
in turn increases the flux of extracted neutrons. 

Neutron monochromator 

Only a narrow range of wavelengths around 8.9 A of incident neutrons is suitable for 
the dominant 1-phonon downscattering in superfluid helium. This allows to benefit from 
a crystal monochromator which reflects apt neutrons and diminishes background at the 
experiment. 

Crystal monochromators reflect neutrons of the desired wavelengths away from the 
primary beam under the d-spacing dependent Bragg angle: 



The d-spacing of a monochromator for A = 8.9 A has to be d > A/2 = 4.45 A. In 
a perfect single crystal the line width for Bragg reflected neutrons is extremely small, 




(4) 
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Figure 2: Schematic sketch of the source with extraction. 1) Intermediate storage 
volume, 2) semidiffuse extraction channel, 3) 70 K heat screen, 4) 4 K heat screen, 
5) vessel for superfluid helium, 6) converter volume. Neutrons entering the channel 

with E± > mgh ( ) are reflected back, the energy components are redistributed 

due to the rough surfaces of the intermediate volume thus it may pass the channel in 
the next attempt. A neutron with E± < mgh ( — ) can pass the channel immediately. 



Ak/k = 10 4 [13] leading to a very narrow acceptance angle of incident neutrons: 

Ah 

— = cot 9A9. (5) 

k 

The incident beam has a divergence of typically ±2° at 8.9 A due to the m = 2 super- 
mirrors used in the guide. Therefore a "mosaic crystal" [13] will be used. Such a crystal 
can be regarded as a collection of microscopic small perfect crystals with differing angles 
e with respect to the overall crystal orientation. Although the angular distribution is in 
general arbitrary, it resembles a cylindrically symmetric Gaussian distribution: 

W(e) = ^e W (-e 2 /27 1 2 ), (6) 
V27r?7 

where 77, called the mosaicity, is the width. To obtain a high acceptance the mosaicity 
should be in the range of the divergence of the incident beam [13]. A high integral re- 
flectivity of the crystal is demanded to ensure an acceptable 8.9 A flux for the converter. 
This requires a small absorption cross-section. 

There are two materials which match these requirements, mica and graphite intercalated 
compounds (GiC) with alkali metals. Due to the small mosaicity « 0.3° of mica we are 
employing a potassium intercalated graphite monochromator of the type previously devel- 
oped at the NIST [15]. The d-spacing for highly oriented pyrolytic graphie is d = 3.35 A 
with a typical mosaicity of 1° — 2°. It is increased by placing guest species in between 
graphite layers (see fig. [3]), this process is called intercalation. 

The monochromator for GRANIT will be made of stage-2 (C24K) potassium inter- 
calated graphite providing a lattice spacing of d = 8.74 A giving a take-off angle of 
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Figure 3: Crystal structure of graphite and staging of GiC. The intercalant diffuses 
into the graphite and thus increases the lattice spacing. Different stages of GiC can 
be produced. The stage number refers to the number of unperturbed graphite layer 
between two layers of intercalant atoms. 



29 = 61.2°. The typical mosaic spread of the produced crystals is 77 = 1.5° — 2.2° 
which matches fine with the incident divergence of beam H172. An integral reflectiv- 
ity of r > 80 % can be achieved [15J. Furthermore the absorption cross section of graphite 
(0.0035 barn) and potassium (2.1 barn) are rather small. The d-spacing of Rubidium and 
Caesium intercalated compounds are of the same magnitude. Rubidium would be an in- 
teresting alternative due to its small cross section of 0.38 barn, whereas Cesium (29 barn) 
is less interesting. For these two alkali metals the production method is not yet sufficiently 
refined. 

Alkali intercalated graphite compounds are conveniently produced using the "two- 
bulb" technique, where the graphite is maintained at a temperature T c which is higher 
than T a of the alkali metal [16J. The stage which is formed depends on the temperature 
difference AT = T c — T a and the quantity of potassium available. For stage- 1, stage-2 we 
are employing 5g, lg ampoules of potassium at T a = 255 ± 3 °C and A7\ = 10 ± 3 °C, 
AT2 = 102 ±3 °C, respectively. A photograph of such a two-bulb cell is shown in figurelH 

A second monochromator with a take-off angle of 29 — 112.5° will include a set of 
stage-1 (CsK) crystals {d = 5.35 A) placed close to the first on a rotary table. The two 
different monochromators can be interchanged (see fig-EJ) making two separate 8.9 A beam 
ports at H172 available. The second beam will feed a position with neutrons for further 
tests and developments on liquid helium based UCN-sources, and later the cryo-EDM 
experiment [9J. 



Converter and cryostat 

As the incident neutron beam has a divergence of a « 2°, the differential flux will be a 
function of position inside the volume. In order to guide divergent cold neutrons incident 
on the walls it is convenient to use polished wall materials with a high critical angle. 
For the storage of UCN we need a high Fermi potential, rough surfaces for an isotrop 
momentum distribution, and a small loss per bounce coefficient [i to minimise wall losses. 
Obviously it is impossible to have polished and rough surfaces at the same place, therefore 
the source will have rough surfaces only on the entrance and exit window of the converter 
volume, whereas the rest is polished. Beryllium, beryllium-oxide, and diamond like car- 
bon (DLC) all have Fermi potentials Vp > 250 neV, and ji < 1 • 10 -4 . 
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Figure 4: Typical cell made from Pyrex used for intercalation. The potassium 
is on the left and causes at T > 250 °C the brown colouring of the glass. On the 
right: a stage-2 intercalated graphite. 




Figure 5: Monochromator setup on rotary table, a) stage-2 potassium intercalated 
graphite monochromator in position, b) stage- 1 (potassium) monochromator. This 
setup allows relatively short switching time. 
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Figure 6: Simulation of UCN-density as a function of converter volume length, 
various rectangular cross sections, and two wall loss coefficients. The incident beam 
is taken to have a Gaussian divergence a > 1.8° dependent on the converging guide 
and differential flux — 3 • 10 8 cm _2 s _1 A _1 at the entrance of the volume. 

These values result from the simulations described in the text. The extraction 
opening is a hole of = 10 mm. UCN once passed through the hole are taken to 
be lost from the volume, in reality there is a certain reflection probability. 



Calculations for a divergent beam, a Fermi potential Vp — 252 neV, a mirror like wall 
(m = 1), and varying loss coefficients are shown in flgure[6l The incident beam flux on the 
monochromator with a reflectivity of r = 85 %, was taken to be d0/dA|A* = 6-10 8 s -1 cm -3 . 
This value was calculated from known cold source data and a transmission simulation 
for the existing guide. The m = 2 neutron guide between monochromator and source 
converges from a 80 x 80mm 2 to the cross section indicated in the figure over a length of 
3 m. Both guide parts have been simulated with a Monte Carlo algorithm. 

The combination of simulated transmissions and divergence for the beam and analytic 
calculations for a converter volume of length 600 mm, a cross section of 70x70 mm 2 , and 
an average \i = 1 • 10 -4 give a density of p calc > 1250 cm 3 . 

The behaviour of the coatings during cooldown still has to be investigated. Thermal 
tension could create cracks and degrade the surface quality significantly. A cooling test 
with several coatings and substrates will be performed to find the best choice. For filling 
and cooling the converter volume with liquid helium we employ a cryostat previously 
developed at Munich jTUl HZ] - Primary cooling power is provided by a Gifford McMahon 
cryocooler with 1.5 W at 4.2 K. It liquefies helium and provides the cooling power for 
the 50 K and 4 K heat screens. The liquefaction is described in detail in ref. [17]. A 
4 He evaporation stage cools the liquid helium below the A-transition point T\ — 2.177 K, 
which then allows us to use a superleak to remove 3 He. The cooling of the converter 
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volume is achieved with a 3 He closed cycle evaporation stage. Using a roots blower pump 
with 500 m 3 /h nominal pumping speed backed by a 40 m 3 /h multiroots pump we were so 
far able to cool the filled converter to 0.7 K. Further improvements, including a new heat 
exchanger between 3 He and 4 He will allow us to reach temperatures below 0.5 K. More 
details about the cryostat will be published elsewhere. 

UCN selection with semidiffuse channel 

The converter volume will be connected to the intermediate volume by a highly polished, 
beryllium, nickel, or DLC coated guide. This additional volume will be made of DLC 
coated aluminium plates. Having rough surfaces it distributes the UCN evenly on the en- 
trance window of the semidiffuse channel and avoids cryogenic complications which would 
arise form a direct attachement of the channel to the converter volume. The optimum 
cross section was determined by simulations with Geant4UCN [18J to be 40x40 mm 2 , the 
length 300 mm is given by the dimensions of the spectrometer. A challenging task is 
to minimise heat input without reflecting UCN back into the converter volume before 
they have reached the intermediate volume. Using thin aluminium or mylar foils as heat 
screens seems to be an option, although losses due to multiple passages will increase due 
to multiple passages. The UCN will be extracted via a 0.2 mm high horizontal semid- 
iffuse extraction channel (see fig. ED. Measurements [12] have shown that such a channel 
increases the storage time and selects only the desired fraction of the phase space. The ex- 
traction channel is made of DLC coated quartz plates [19], the lower surfaces are polished, 
the upper rough. The channel dimensions are h = 200 /im, / = 100 mm, w = 300 mm, 
which will provide a reflectivity r > 85 % for neutrons with a vertical energy component 
E± > mgh = 20 peV. 

Conclusion 

The crystals of the monochromator, the converter, and the semidiffuse channel have all 
been tested separately. Their integration remains a challenging task. Present calculation 
with an incident differential flux d(f)/d\\\* = 6 • 10 8 s -1 cm -3 on the monochromator give 
a UCN density of Pucn ~ 1250 cm -3 in the converter, p int = 800 cm -3 in the intermediate 
volume. This yields for a critical velocity of 7 m/s for the materials used an available 
phase-space density of r He ~ 0.18 cm -3 (m/s) -3 . Compared to the phase-space density 
r Tu rbine ~ 0.013 cm _3 (m/s) -3 of the UCN turbine at the ILL this is more than a factor 
ten of improvement. This calculation assumes perfect conditions and in all parts optimal 
transmissions, experience with other simulations have shown that the actual value might 
be lower. 
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